Doping front instabilities in organic semiconductors: a means for optimizing 

optoelectronic devices 
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Recently, it was demonstrated that electrochemical doping fronts in organic semiconductors ex- 
hibit a new fundamental instability growing from multidimensional perturbations [Phys. Rev. Lett. 
107 016103 (2011)]. In the instability development, linear growth of tiny perturbations goes over 
into a nonlinear stage of strongly distorted doping fronts. Here we develop the nonlinear theory 
of the doping front instability and predict the key parameters of a distorted doping front, such as 
its velocity, in close agreement with the experimental data. We show that the instability makes 
the electrochemical doping process considerably faster. We obtain the self-similar properties of the 
front shape corresponding to the maximal propagation velocity, which allows for a wide range of 
controlling the doping process in the experiments. The theory developed provides the guide for 
optimizing the performance of organic optoelectronic devices. 



I. INTRODUCTION 



The discovery and development of 7r-conjugated poly- 
mers - plastics with properties of semiconductors, often 
termed organic semiconductors (OSCs) - has revolution- 
ized the field of electronics, and was awarded with the 
Nobel Prize in Chemistry in 20001"^. The OSCs, in con- 
trast to crystalline inorganic semiconductors, are sim- 
ple to process and relatively inexpensive in production, 
while their soft and conformable character enables large 
number of novel technological applications ranging from 
optoelectronics and photovoltaics to sensors and artifi- 
cial muscles^"—. The unique combination of electronic 
and ionic conductivity in OSCs introduces the exciting 
functionality that triggers experimental and technologi- 
cal activities in the nascent field of so-called organic (i.e. 
based on OSCs) electronics. At the same time, as for- 
mulated in Ref ,5, " the understanding of the fundamental 
processes that take place [in OSCs] is still in its infancy" , 
especially if compared to the theoretical fundament of the 
field of crystalline inorganic semiconductors-'^'^'U-. The in- 
adequate understanding of fundamental OSC properties 
limits the technologies to " cut-and-try" methods and de- 
mands development of the theory on the subject. In line 
with this demand, the goal of the present work is to attain 
better fundamental understanding of the doping process 
in OSCs, which may provide a guide to optimizing per- 
formance of optoelectronic devices such as, for example, 
light-emitting electrochemical cells. 

In the process of electrochemical doping the OSCs 
can be transformed from an almost insulating state to a 
metallic-like state with the conductivity increased by sev- 
eral orders of magnitude^ ' ^^d^ . Apart from conductivity, 
this doping transformation also allows for tuning other 
key properties of OSCs, such as color, volume, surface 
energy and photoluminescence capability, in a control- 
lable way. This possibility can be profitably employed in 
immerous applications, e.g., in optoelectronics"^^^. Oper- 
ationally, the electrochemical doping can be performed 
in-situ by applying a potential difference between two 
metal electrodes in contact with an OSC, blended with 



a solid-state electrolyte, as illustrated in Fig. 1 (a). Be- 
cause of the applied electric field, positive and negative 
ions redistribute in the active material and allow the in- 
jection of electrons (n-doping) and holes (p-doping) into 
the OSC from the cathode and anode, respectivelyi^"— . 
As the electronic charges populate the OSC, its conduc- 
tivity increases by several orders of magnitude in com- 
parison with the initial conductivity supported by ions. 
It has been demonstrated experimentally that the dop- 
ing process develops in the form of two counter propa- 
gating fronts illustrated in Fig. 1 (a), which separate 
the regions of p- and n-doped material from the undoped 
onei2r2S. When the two fronts meet, a p-n junction is 
formed in which the recombination of holes and electrons 
results in the visible light^iil"^. Thus, the velocity of 
doping front propagation determines the turn-on time of 
the light-emitting electrochemical cells. 

As pointed out above, the present understanding of 
the doping process in OSCs is far from being ade- 
quate, with only few theoretical papers focusing on the 
subjec t^^i^^" — . In particular, the analytical formula for 
the velocity of a planar doping front has been found only 
recently in Refi^, which demonstrated good agreement 
with the experimental observations. The analytical the- 
ory has shown that the doping rate is controlled by the 
mobility of heavy ions that results in a slow turn-on time 
of light-emitting electrochemical cells. At present, one 
of the most difficult problems to solve is to speed-up 
the doping process and to reduce the turn-on time. The 
speed-up could be achieved with the help of a new ba- 
sic phenomenon - the multidimensional doping front in- 
stability - discovered recently in OSCs in Ref^iS. This 
new fundamental instability distorts the shape of dop- 
ing fronts, increases locally the electric field ahead of 
the fronts and hence increases their propagation veloc- 
ity. The recent paper^- has developed a linear theory 
and provided the experimental evidence of the instabil- 
ity; it has also proposed a possible way to amplify the 
instability by means of minor humps introduced at the 
electrode surfaces. Figure 1 (b) illustrates the strongly 
corrugated doping fronts obtained in the experiments of 
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and plasmas^""—. In particular, in the present work we 
discuss similarity and difference between the doping front 
instability in OSCs and the hydrodynamic instability of 
a flame, i.e. a slow combustion front^^"^^. The electro- 
chemical doping process provides a fascinating opportu- 
nity to observe directly and investigate a number of inter- 
esting nonlinear phenomena: front propagation, multi- 
dimensional instabilities and pattern formation, with sur- 
prisingly different outcome for p- and n-fronts. In addi- 
tion, the present analysis suggests different regimes of 
the doping front instability for the 2D geometry of OSC 
films and for the 3D geometry of the bulk optoelectronic 
devices^: thus, the nonlinear dynamics of the instability 
depends also on the system dimension. 

The paper is organized as follows: In Sec. II we remind 
basic results on the linear stage of the doping front in- 
stability needed for the subsequent nonlinear analysis. In 
Sec. Ill we derive the nonlinear equation for the p-doping 
front bent by the instability. In Section IV we analyze 
stationary periodic solutions to the nonlinear equation. 
In Sec. V we compare the theoretical results to the ex- 
perimental data and discuss technical implications of the 
theoretical predictions. The paper is concluded by a brief 
summary in Sec. VI. 



II. LINEAR STAGE OF THE DOPING FRONT 
INSTABILITY 



FIG. 1: Schematic of a light-emitting electrochemical cell (a). 
Experimental photographs of Ref.^^ showing the time evolu- 
tion of the doping fronts in a light-emitting electrochemical 
cell with the corrugations imposed at the anode surface (b). 
The corresponding time instances are indicated on the pho- 
tographs. 



Refji^. Still, the dynamics of the strong instability at 
the nonlinear stage remained unclear and demanded deep 
theoretical studies. 

In the present paper we develop the nonlinear theory of 
the doping front instability, focusing on the velocity and 
shape of the corrugated fronts. For the 2D geometry of 
the OSC films, the theory predicts increase of the doping 
front velocity by a factor of two in comparison with the 
planar front velocity, which is in a good agreement with 
the previous experimental data. We obtain self-similar 
properties of the front shape corresponding to the maxi- 
mal velocity, which promises a wide range of controlling 
the doping process in the devices, e.g., through the de- 
sign of the electrodes. We show that moderately corru- 
gated doping fronts are preferential for achieving higher 
velocity as compared with the strong corrugations ob- 
tained in Rcf.^^ and illustrated in Fig. 1. The developed 
theory provides a guide for optimization of organic opto- 
electronic devices utilizing electrochemical doping. 

In addition, it is expected that doping fronts in OSCs 
will become a new popular object of nonlinear science 
alongside with other fronts and interfaces in fluids, gases 



In this section we remind basic results of Ref-i^^ on the 
doping front instability at the linear stage, when pertur- 
bation amplitude is infinitesimal. Details of the linear 
analysis are needed to understand the nonlinear theory, 
which we develop below in Sec. Ill and IV. The instabil- 
ity arises because speed of a doping front depends on the 
electrical field ahead of the front as it was derived in2^. 
We present the p/n-doping front velocity in the vector 
form as 



Up,„ = ±^(/i++^_)E. 

nh.e 



(1) 



Here E ~ — V(/> is the electric field intensity just ahead 
of the front, (f) the electric field potential, no the initial 
ion concentration in the active material, rih^e the final 
concentration of the holes/electrons behind the front de- 
termined by the thermodynamic properties of a partic- 
ular OSC, /i+ and ^_ the mobilities of the positive and 
negative ions, respectively. The "±" sign in Eq. ([T]) in- 
dicates that the p-type {+) and n-type (— ) doping fronts 
propagate in opposite directions. We start our analysis 
with a single p-doping front propagating stationary with 
velocity Up in a constant uniform external electric field 
E — slxEq pointing along x-axis, where a.^ is the unit 
vector in x-direction. The analysis of the n-front insta- 
bility may be performed in a similar way. In the case 
of a planar stationary front, the general 2D equation for 
the p-front position x = Xp{t, y) reduces to Xp[t) = Upt. 
Next, we consider small perturbations of the stationary 
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FIG. 2: Schematic of the doping front instabiUty. 



p-doping front 



Upt 



(2) 



where Xp cx exp{at + iky) stands for one Fourier mode 
at the Unear instabiUty stage, k = 27r/A is the pertur- 
bation wave number characterizing infinitesimal bending 
of the doping front (A is the perturbation wavelength), 
a is the instability growth rate, see Fig. 2. The purpose 
of the linear stability analysis is to find the dispersion 
relation a{k). In the case of Re((T) > valid for a certain 
domain of perturbation wave numbers we conclude that 
the planar front is unstable with respect to multidimen- 
sional perturbations. Small bending of the doping front 
induces also perturbations of the electric field in the un- 
doped region y, x) cx exp(at+iky). The doped region 
may be treated as equipotential due to high conductivity. 
Thickness of the doping front may be characterized by a 
length scale Lp, of about 10"^— 10~'^mm, determined by 
ionic diffusio n^'^1^^ . Since the characteristic size of the 
experimentally observed p-front perturbations (approx- 
imately 0.2mm) is much larger than Lp, we may treat 
the front as infinitesimally thin. Most of the time the 
doping fronts are sufficiently far away from each other, 
k{Xn — Xp) ^ 1, therefore, the instability of one front is 
not affected by the other. The linearized form of Eq. ([T]) 
may be written in components as 



dXp 
dt 



no 



(M4 



dx ' 



no ^ 



ikUpXp = — {fi+ + fi^) — . 

nh oy 



(3) 



(4) 



Solving the Laplace equation = for perturbations 
of the electric potential in the udoped region 



dx"^ 



kU = 0, 



we find the solution 

(j) cx exp((Tt -I- iky — \k\x). 



(5) 



(6) 



Eqs. dS]), dH), ([6]) may be reduced to a single linear equa- 
tion describing the perturbation growth with time 



^ = Up\k\Xp. 



(7) 



Thus, we obtain an instability for the p-type doping front, 
where the initially small perturbations grow exponen- 
tially as Xp cx exp((T<), with the growth rate 



;7p|fc|. 



(8) 



The same result holds for the n-type doping when Up 
is replaced with [/„. The obtained dispersion relation 
cr cx I A: I is mathematically similar to the Darrieus-Landau 
(DL) instability encountered in combustion, astrophysics 
and laser fusion22ri^. The role of the electric field in the 
new instability may be also compared to that of the gas 
velocity field in flames. 

According to the above analysis, the doping front in- 
stability develops because of a positive feedback of two 
effects: First, the local doping speed depends on the elec- 
tric field, Eq. dl]); and, second, the electric field increases 
close to the humps at a curved surface of the conducting 
materials^,. The latter is reminiscent of the St. Elmo's 
fire phenomenon, which visualizes strong increase of elec- 
tric field at sharp conducting surfaces by numerous dis- 
charges in the surrounding air. In the present configu- 
ration of a curved doping front, any small hump at the 
front produces local increase of the electric field in the 
undoped material close to the hump, which, in turn, fa- 
cilitates further growth of the hump, thus leading to the 
instability. 



III. THE NONLINEAR EQUATION FOR THE 
P-DOPING FRONT BENT BY THE 
INSTABILITY 

A smooth shape of the p-doping front observed in the 
experiments^^ allows reducing the full diffusion-mobility 
model^^- to a single nonlinear equation. Similar to the 
previous section, we consider a single p-doping front 
propagating in the initially uniform external electric field 
E — SlxEq. In developing the nonlinear theory of the 
doping front instability for perturbations of considerable 
amplitudes, we have to take into account interaction of 
modes with different wave numbers. Then Eq. ([7]) has 
to be presented in the form 



dXp _ 

-gf - UpJXp, 



(9) 



where the DL operator J implies product by \k\ in Fourier 
space, i.e. 



JXp{t,y)= / \k\Xpk{t,k)eM^ky)dk (10) 
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in the 2D geometry of an OSC film. Equation ([9|) will 
be obtained below as the principal linear terms in the 
nonlinear equation for the curved p-doping front. 

We study dynamics of the p-doping front as a geomet- 
rical surface 



Gir,t)=x-Upt-Xp{t,y) 
using the identity 







dr 



VG- 



dG 



0, 



dG 

'dt dt ' " ' dt 
which leads to the normal front velocit v^"'^^ 



dr 

'dt 



1 dG 
\WG\~dt 



Up + dtXp 
1 + idyXp) 



(11) 



(12) 



(13) 



Here we have introduced the unit normal vector to the 
surface 



|VG| 



1 -I- idyXpY 



(14) 



On the other hand, the normal velocity of the curved p- 
doping front is determined by the electric field just ahead 
of the front, see Eq. ([1]), which leads to 



Up + dtXp 



no 



1 



{dyXpY "'^ 



-^{^i+ + ^l.){k^Eo-Vct>). (15) 



Taking into account the relation between the planar front 
velocity Up and the external electric field Eq, we write Eq. 
psp in components as 



U, 



±^^Up-^d., 

{dyXpY ' EO 



dyXp 



u„ 



1 + {dyXp 



dtXp _ Up ■ 



(16) 



(17) 



We also employ the solution to the Laplace equation for 
one Fourier mode of the perturbed electric field cf) cx 
exTp{iky ~ \k\x) and find the general identity 



dj^^ — — J(j) = —J{dy(j)/ik) = —Hdyq 



(18) 



where the Hilbert operator H means multiplication by 
\k\/ik in Fourier space. Then Eqs. (IT6)) . ([T7| may be 
reduced to 



Up^dtXp 



{dyXp 



1 + {dyXpf 



= H 



dyXp 



~ l + Up^dtXp 

{dyXpY 



1 



(19) 



Within the limit of weak second-order nonlinearity we 
simplify Eq. as 



u: 



.i dXp 
dt 



JX„ 



fdXp^ 


2 

) +Up'H 


'dXp 


dXp 






dt 


dy 



(20) 




FIG. 3: Schematic of the Huygens nonlinear stabilization: 
solid lines represent successive front positions; dotted arrows 
indicate local direction of front propagation. 



The first two terms in Eq. (|20p describe the linear in- 
stability growth similar to Eq. (jH). The third term cor- 
responds to the nonlinear Huygens stabilization of the 
instability due to the front propagation, which is well- 
known, e.g., in combustion science^Sr^. Because of the 
Huygens stabilization convex parts of the front (humps) 
become smoother, while concave parts converge to sharp 
cusps, as illustrated in Fig. 3. Propagation speed of the 
cusps increases in comparison with the planar front ve- 
locity Up, and the cusps tend to catch up with the humps. 
When linear growth of the humps because of the insta- 
bility is balanced by the nonlinear increase of the cusp 
velocity, we should expect formation of a curved station- 
ary or quasi-stationary p-doping front propagating with 
an increased velocity in comparison with Up. Still, one 
should employ the similarity between combustion fronts 
and the doping fronts in OSCs with caution, since these 
two phenomena demonstrate also a number of important 
different features apart from similar ones. In particular, 
the normal velocity of a thin combustion front may be 
typically treated as a constant value independent of the 
large-scale gas fiow of the fuel mixture22r^. On the con- 
trary, the normal velocity of a doping front is determined 
by the electric field in the undoped material just ahead of 
the front, see Eq. ([1]). Finally, the last term in Eq. (^0)) 
has not been encountered before in combustion science or 
any other study of front dynamics. This term is intrinsic 
for the doping fronts in OSCs, and, as we will see below, it 
provides additional nonlinear destabilization of the front 
since it has the same origin as the linear destabilizing 
term JXp. The nonlinear destabilization may be quite 
strong in the 3D geometry so that the whole approach of 
weak nonlinearity fails and one has to go back to the fully 
nonlinear equation Eq. (|19p . which is much more difficult 
for solution than Eq. Still, in the 2D geometry of 

the OSCs films studied experimentally in Ref.^^ the ap- 
proach of weak nonlinearity provides reasonable results 
able to explain the experimental data both qualitatively 
and quantitatively. 
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In the studies of the unstable p-doping fronts we also 
have to take into account linear stabilization of the in- 
stability at short wavelengths below the cut-off Xp oc Lp. 
The short wavelength stabilization is related to the small 
front thickness and, therefore, was not accounted in 
Eq. ([T]). Though negligible at the experimentally ob- 
served length scales, the stabilization prevents the ul- 
traviolet divergence of the perturbation spectrum in the 
doping front dynamics by use of the modified dispersion 
relation 



a^U„ \k 



27r 



(21) 



The above dispersion relation may be also presented in 
the form of a linear integral-differential equation similar 
to Eq. (HD 



^i dXp 
dt 



Xp d'^Xp 
2^ dy^ ■ 



(22) 



The linear stabilization terms have the physical mean- 
ing of the local front velocity depending on the front 
curvature22r^; these terms may be also incorporated as 
curvature corrections to the doping front velocity in Eq. 
([T]). Then, taking into account the linear stabilization 
term in the nonlinear equation Eq. (j20p . we come to the 
final equation for the p-doping front 



1 dXp 



J x„ 



dXp 
dy 



1 



-H 



dXp dXp 



dt dy 



+ 



Xp d^Xp 



2tt 



Equation has some 



dy^ 
(23) 

common features with the 
Sivashinsky equation and other nonlinear equations de- 
rived within the theory of the DL hydrodynamic flame 
instability^i^i"— . Due to this similarity, we may em- 
ploy the mathematical methods of investigating the DL 
instability to the problem of doping fronts in OSCs as 
demonstrated in the next section. 



IV. SHAPE AND VELOCITY OF STATIONARY 
CURVED P-DOPING FRONTS 

Here we consider stationary periodic solutions to 
the derived nonlinear equation Eq. (j23p in the form 
Xp{y,t) — tAU + f{y), where AU is additional increase 
of the doping front velocity because of the instability 
and f{y) is a periodic function of the period A. Despite 
the imposed limitations, the stationary periodic solutions 
provide the clue to the general understanding of the p- 
doping front dynamics. We reduce Eq. (^5]) to the sta- 
tionary form 



AU/Up = (1 + AU/Up) Jf + 



dl 
dy 



27r 9y2 ■ 



(24) 



We point out that Jf is the key term describing the 
front instability at the linear stage; thus the product 



{AU/Up) jf with AU > means additional nonlinear 
destabilization of the front. It should be stressed that 
such a term has not been encountered before in the com- 
bustion studies of the DL instability, and it modifies the 
nonlinear dynamics of the p-doping fronts quite strongly 
as compared to the combustion fronts. 

By employing the method of pole decomposition^"—, 
we look for A-periodic solutions to Eq. in the form 



A 



2N 

p Y^i 



2 iy-va) 



(25) 



where the complex values ya represent the poles, un- 
known so far. Substituting the pole solution Eq. ((25)) 
to the stationary equation (1^^ we derive the set of equa- 
tions for the poles ya- In these calculations, action of the 
DL operator J on Eq. ([^5]) requires additional comments. 
The definition of the operator J may be re- written as fol- 
lows 

J if (y)] - ^ E ^'^ 7 '^os [^^'^ (y - 2/')]/ (y) dy', (26) 



where the integral in Eq. ([26| is similar to the Fourier 
transform. We integrate Eq. ([25)) by parts using the 
following identities^ 



smy 



cosh z — cos y 



2E 



e cos ny 



sinh; 



n=l 



cosh z — cos y 



'sinny, z > 0, (27) 



1, z > 0. (28) 



As a result, we can specify action of the DL operator on 
the pole solution Eq. ([25)) as 



Jf = k 



k sinh [fcealm(j/a)] 



cosh [kealvaiya)] - cos [k (y -I- lm{ya))] ' 



(29) 



with the designation Sa = sign [Im(?;ct)] introduced for 
brevity. Then, we substitute Eq. to Eq. ([M]) and 

perform straightforward calculations to find the set of 
equations for the poles 



Xp/X 



2N-1 

; -, >^ cot 

1 + AU/Up 



ya - yp 



X 



isign(Imya) = 0, 



(30) 

where N pairs of poles are involved in the solution. In 
the derivation of Eq. ([30)) we also employed the following 
trigonometric formula 



cot a cot b = cot (a — 5) [cot b — cot a] — 1 . 



(31) 



Apart from the equations for the poles, the calculations 
yield an equation for the velocity increase 



AU/Up 
1 + AU/Up 



A 



1 



NXp/X 
1 + AU/Up 



(32) 
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FIG. 4: Velocity increase AU/Up of the p-doping front for 
stationary periodic solutions versus the dimensionless period 
A/Ap: (a) for one pair of poles N = 1; the solid line indicates 
the limited solution, the dashed line indicates the diverging 
solution; (b) for different numbers of pole pairs A*'; the solid 
lines indicate maximal possible velocity increase for a fixed 
period, the dashed lines present other solutions. 

where N is an integer standing for the number of the pole 
pairs. Solution to Eq. (15^ is illustrated in Fig. 4 (a) for 
iV = 1: the solid lines present the solutions with velocity 
increase limited for any Ap/A; the dashed lines show the 
solutions with unlimited velocity increase. The latter, 
presumably, does not have a physical meaning and hence 
will be omitted in the following. The solution presented 
by the solid line is determined by a simple formula 

AU/Up = NXp/X for N < A/Ap, (33) 

which follows from Eq. ([5^. We can also obtain the 
same result for the velocity increase by taking the spatial 
averaging of Eq. (IM)) as 

AU 1 } fdfV , , , 



The velocity increase is presented in Fig. 4 (b) versus 
the scaled period of the solution A/Ap for different num- 
bers of the pole pairs. As we can see from Fig. 4 (b), 
a curved p-doping front is possible only for a sufficiently 
large period A exceeding the instability cut-off A > Ap; 
the branch starting at A = Ap corresponds to one pair of 



poles, = 1. The solution with N pairs of poles becomes 
possible for A > iVAp and it provides the maximal front 
velocity within the interval of iV < A/Ap < + 1. As 
we increase the period A in comparison with Ap, multiple 
stationary periodic solutions of the same period become 
possible with different numbers of the pole pairs. In Fig. 
4 (b), the solid lines indicate the solutions correspond- 
ing to the maximal velocity increase for a fixed period A; 
other solutions are presented by the dashed lines. The 
solutions with maximal velocity involve also the largest 
possible number of poles for a fixed period. Presum- 
ably, the solutions with the maximal propagation velocity 
and the maximal number of poles are stable, while other 
solutions (shown by the dashed lines) are unstable. In 
the context of combustion fronts, such a result has been 
demonstrated in Refi^ by investigating stability of the 
pole solutions; we expect the same conclusion to hold for 
the doping fronts in OSCs as well, though the detailed 
stability analysis requires much work and is left for the 
future. Figure 4 (b) demonstrates also that the maximal 
velocity increase tends to the limiting value 

maxAC//C/p = 1 (35) 

as the period A becomes much larger than the cut-off 
Ap, i.e. for A/Ap — oo. The limit of large scales in 
comparison with the cut-off is of the most importance 
for the analysis of the experimental results^^. Thus, we 
expect strongly curved p-doping fronts in OSC films to 
propagate twice as fast in comparison with the planar 
fronts. 

Next, by using Eq. we investigate the front shape 
predicted by the stationary solutions. Translational in- 
variance of the solution Eq. (pS)) implies that we can 
chose arbitrary Re(?/Q); hence we set Re(?;c<) = 0, so that 
the poles ya become pure imaginary = ±iBa. Tak- 
ing into account complex conjugate pairs, we rewrite Eq. 
(|25p in a more convenient form 

A ^ 1 

/ (y) = ^ E In 2 ^^^^^ " ^^y^^ ^^^^ 

and the set of equations for the poles, Eq. ([5(1)) . reduces 
to 



^ (^1 + = coth [kBn] + (37) 




for n ~ 1,...N. Figure 5 (a) shows positions for the 
poles Ba calculated for Eq. (|37p for different number of 
pole pairs A^ and the solution period A/Ap = A^ corre- 
sponding to the maximal velocity increase. As we can 
see in Fig. 5 (a), scaled positions of the poles B^ de- 
crease strongly for any fixed a as the number of poles 
increases. Since the solution (I36p involves linear super- 
position of the poles, then we can study relative roles of 



FIG. 5: Illustration of the pole solutions: (a) Scaled pole po- 
sitions Be versus the number of pole pairs for A'^ = X/Xp; (b) 
Solution for three pole pairs for X/Xp = N = 3 and separate 
contributions of different pairs, q = 1, 2, 3 selected in the plot 
(a). 



different pole pairs separately. As an illustration, Fig. 5 
(b) presents the solution with three pairs of poles 
for A/Ap = TV = 3 together with separate contributions 
from different poles a = 1 — 3 (marked in Fig. 5 (a)). We 
observe in the figure that the pole pair with the minimal 
number a = 1 plays the dominant role in building the 
amplitude of the final hump, while the poles with higher 
numbers are responsible mostly for fine details. In that 
sense, the role of poles in the solution may be com- 
pared qualitatively to the role of harmonics in Fourier 
series. 

Substituting the solutions for the poles into Eq. (j36p 
we find the shapes of the curved p-doping fronts. In 
agreement with the experimental observations. Fig. 1, 
the shapes are built of smooth humps facing the undoped 
material and sharp cusps pointing to the doped one, e.g. 
see Fig. 6 (a) with the solution plotted for A/Ap = 2. 
For illustrative purposes. Fig. 6 (a) presents three peri- 
ods of the solution. We remind that in this section we 
focus on the stationary solutions with humps and cusps 
reproduced periodically. Still, these humps provide the 
key to understanding irregular front patterns as well like 
those observed experimentally and shown in Fig. 1. Fig- 
ure 6 (b) presents the stationary humps corresponding 
to the maximal front velocity and the maximal number 
of poles found for different values of the scaled period 
A/Ap = = 1 — 10. In Fig. 6 (b) we observe only mi- 
nor variations of the hump shape with increasing period. 



FIG. 6: Shapes of the curved p-doping fronts: (a) periodic 
solution plotted for A/Ap = 2. (b) Shapes of humps for 
A/Ap = iV = 1,2,5, 10. 



The variations concern mostly the cusps, which become 
sharper as we increase the hump size A in comparison 
with the instability cut-off Ap. We remind that the limit 
of large length scales, A/Ap — oo, is the most interesting 
from the experimental point of view. As we can see in 
Fig. 6 (b), the shape of a hump becomes self-similar in 
the limit of large length scales, and this shape may be ap- 
proximated with a very good accuracy by a solution with 
a limited number of pole pairs. In particular. Fig. 6 (b) 
demonstrates negligible difference between the solutions 
obtained for X/ Xp = N — 5 and 10. 

At the end of this section we point out that the pole 
solutions may be written in an explicit analytical form 
for small numbers of pole pairs A^ = 1,2. For one pair of 
poles (A = 1) we find from Eq. §7} 



A 

A„ 



AU \ 



coth [kBn 



which leads to the front shape 



27r 



1 - Wl - 



Ap A / 27ry 
- — ^ cos — — 

X + Xp) \ X 



(38) 



(39) 



In the same way, for two pairs of poles (A^ = 2) we obtain 

3Ap/A 




(40) 
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and the front shape 



■' 2-K 4 



l + A/2 



1 - A 



(41) 



/-V2 - A /27ry\ 

V2— — cos — — +COS , , , 

1-A \ ^ ) \ ^ ) 



2 



where the designation A 
introduced for brevity. 



[3Ap/(A + 2Ap)]" has been 



V. COMPARISON TO THE EXPERIMENTS 
AND TECHNICAL IMPLICATIONS 

In this Section we compare the theoretical predictions 
to the experimental results of Refr^^ and demonstrate 
how theoretical understanding may be employed for bet- 
ter control and optimization of the electrochemical dop- 
ing process in optoelectronic devices. In particular, Fig- 
ure 7 compares positions of the doping fronts as predicted 
by the theory and measured experimentally. Here, the 
solid lines present front positions assuming the planar 
front shape (labeled "theory 1") and accounting for the 
instability and the curved shape of the p-doping front (la- 
beled "theory 2"). The markers show average positions 
of the strongly corrugated fronts obtained experimentally 
in Ref.^^ and illustrated in Fig. 1; the error-bars indicate 
average spreading of the front brush. Details of the theo- 
retical calculations for the geometry of the planar fronts 
may be found in Rcf.^^. Here we remind two important 
features of the process taken into account in the calcu- 
lations: (i) according to the experiments^, the n-doping 
front starts later and propagates noticeably slower than 
the p-doping front; (ii) the fronts accelerate towards each 
other because of much higher conductivity of the doped 
material in comparison with the undoped one. Due to 
the relatively high conductivity of the doped material, 
the potential difference applied to the electrodes drops 
effectively at the gap between the fronts. As the distance 
between the p- and n-fronts decreases, the electric field 
in the gap increases, which leads to increase of the front 
velocities according to Eq. ([T]) . Theoretically, the electric 
field and the planar front velocities diverge as the fronts 
meet; the front propagation comes to a complete halt at 
that point. The theory of Sec. IV predicts that, in the 
experimental conditions, the doping front instability and 
the corrugated shape may increase the propagation ve- 
locity of the p-doping front by a factor of 2. We employed 
this theoretical result for calculating the average position 
of the corrugated p-doping front labeled as " theory 2" in 
Fig. 7. Since the theory is valid only for p-doping fronts, 
then we use the model of a planar n-doping front as in 
the previous calculations. Still, the n-doping front sweeps 
only about 10% of the active material in the described 
experiments, and, therefore, a particular model for the n- 
front dynamics produces only slight modifications of the 
results. 




Time (s) 

FIG. 7: Position of the doping fronts versus time as predicted 
by the theory for the planar fronts (theory 1) and taking into 
account the instability (theory 2). The experimental results 
of Ref.i^ for the average front positions are presented by the 
markers. The error bars indicate the average spreading of the 
front brush. 



Figure 7 shows that the well-developed front instability 
and the curved shape of the p-doping front may speed- 
up the doping process considerably thus improving per- 
formance of the optoelectronic devices. The theoretical 
predictions for the p-front dynamics are in a good agree- 
ment with the experimental data during most part of 
the doping process (for the time instants up to ^ 65 s). 
The front velocity found experimentally is only a little 
lower than the theoretical expectations for the maximal 
front velocity within that interval. Some deviations in 
the front position and velocity between the theory and 
the experiments observed in Fig. 7 are quite natural, 
since the periodic solutions studied theoretically in Sec. 
IV are not identical to the irregular pattern of Fig. 1 
obtained experimentally. As we discuss below, the differ- 
ence in the front shape should lead to larger velocity of 
the periodic curved fronts of Sec. IV in comparison with 
the front pattern of Fig. 1. The deviations between the 
theory and the experiments become considerable at the 
end of the process, for the time instants about (70 — 80) s. 
On the other hand, the very approach of the average front 
position is misleading at the end of the doping process, 
since the gap between the p- and n-fronts becomes com- 
parable to the size of the humps. 

Figure 8 compares shape of the humps at the curved 
p-doping front as predicted by the theory for a stationary 
periodic solution and obtained in the experimentsA^. In 
order to perform the comparison, we shift and re-scale the 
humps labeled in Fig. 8 (a) keeping equal scales in both 
directions. We point out that comparison is possible only 
for the top part of the humps, since the bottom parts 
are strongly affected by particular irregular conditions 
of front initiation. Figure 8 demonstrates a reasonable 
agreement between the theory and the experiments. Still, 
we can see that the periodic stationary humps studied in 
Sec. IV are somewhat smoother and the cusps in between 
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FIG. 8: Shape of the humps obtained in the experiments^^ for 
the time instant of 54 s, (a), and comparison of the experimen- 
tal results for different humps to the theoretical predictions for 
a stationary periodic solution, (b). The experimental results 
in the plot (b) are shifted and re-scaled in order to demon- 
strate self-similarity of the humps. 



them are more shallow than those observed in the exper- 
iments. In part, the difference may be attributed to the 
limited accuracy of the assumption of weak nonlinearity 
used for the derivation of Eq. from the fully non- 

linear equation Eq. (jl9p . However, the fully nonlinear 
equation Eq. ()19p is much more difficult for the analysis 
than Eq. (|23|) : for this reason, investigation of Eq. (fT9| is 
left for future studies. Large amplitude of the humps in 
Fig. 1 may also originate in a particular way of initiating 
the instability in the experiments-''?. As demonstrated in 
Fig. 1, the strong instability of the p-doping front has 
been initiated by introducing minor corrugations at the 
anode surface. Because of the local increase of the elec- 
tric field at the corrugations, the doping process started 
faster at the corrugation tips. For a while, the regions 
of p-doping developed separately in the form of isolated 
humps close to the anode corrugations without forming 
a single continuous p-front. When the front was finally 
formed, the employed way of initiating the instability 
produced extremely deep cusps - chasms - in between the 
humps, which remained quite deep till the end of the dop- 
ing process thus reducing the average front position and 
the average front speed. Here we would like to stress an 
important difference between the doping fronts in OSCs 
and combustion fronts one more time; the similarity be- 
tween these two phenomena is limited and may be often 
misleading. In the case of combustion, normal velocity of 
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FIG. 9; Absolute value of the electric field (a) and potential 
(b) calculated in the undoped region for the corrugated p- 
doping front of Fig. 8 (a). 
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FIG. 10: Absolute value of the electric field (a) and poten- 
tial (b) calculated in the undoped region for the moderately 
corrugated p-doping front obtained in Sec. IV with the same 
number of humps as in Fig. 9. 



the front propagation is approximately constant, so that 
deep chasms at a combustion front imply strong increase 
of the front speed, and the Huygens stabilization mech- 
anism leads to fast leveling of the chasms as shown in 
Fig. 3. In contrast to combustion, such effects do not 
happen at the doping fronts in OSCs. The doping front 
velocity is determined by the electric field just ahead of 
the front, which is strongly reduced in the chasms. In 
particular. Fig. 9 presents the electric field (a) and po- 
tential (b) calculated for the corrugated doping front of 
Fig. 8 (a). The electric field in Fig. 9 varies by 3 — 4 
orders of magnitude from ^ 12kV/m at the tips of the 
humps to ~ (10 — 100)V/m at the bottom of the chasms. 

For this reason, the Huygens stabilization mechanism 
is much weaker for the doping fronts, and deep chasms 
reduce the average front velocity instead of increasing it. 
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As a result, moderately corrugated doping fronts studied 
in Sec. IV provide higher average propagation velocity 
than the strongly corrugated fronts of Fig. 1. The elec- 
tric field (a) and potential (b) in the undoped material 
ahead of a moderately corrugated doping front studied in 
Sec. IV is shown in Fig. 10; in that case the electric field 
between the humps and the cusps varies much less than 
in Fig. 9, approximately by an order of magnitude, from 
~ 5kV/m to ~ 0.5kV/m. Thus, in order to obtain fur- 
ther increase of the doping rate in comparison with the 
experiments^^, one should avoid extremely deep chasms 
when initiating the doping front instability. This can be 
achieved, for example, by introducing smoother corru- 
gations at the electrodes imitating the curved periodic 
doping fronts obtained in Sec. IV. 

Another important feature of the process is the self- 
similar shape of the humps at the corrugated p-doping 
front, which allows additional control of the front corru- 
gations and, hence, of the final light-emitting p-n junc- 
tion curve. As we have shown in Sec. IV, in the 
limit of large length scales A/Ap ^ 1 inherent to the 
experiments^, the maximal velocity increase tends to 
the limiting value AJ7/C/p — >■ 1 independent of a par- 
ticular length scale A, and the humps attain a self- 
similar shape presented in Fig. 6. We point out that 
the sample sizes in the experiments-^, namely, breadth 
of 1 mm and length of 1 cm, were many orders of 
magnitude larger than characteristic width of the dop- 
ing front, Lp (10~^— 10~'^)mm, determined by ionic 
diffusion^'^'^^. For this reason, humps of both large and 
small size in comparison with the sample sizes may satisfy 
the limit X ^ Xp (x Lp and provide maximal velocity of 
front propagation. As an illustration, Fig. 11 shows two 
p-doping fronts with two different periods of the humps, 
A = 0.3 mm and 1.2 mm, which propagate with the same 
maximal velocity. Visually, front (b) demonstrates strong 
corrugations, while front (a) is almost planar in aver- 
age when compared to the sample sizes. Respectively, in 
the case (b) one should expect a strongly curved light- 
emitting p-n junction, while in the case (a) the p-n junc- 
tion will be almost planar with slight wrinkles, though 
both cases demonstrate the same maximal doping front 
velocity increased because of the instability. 

Finally, we discuss implications of the present analysis 
for the 3D geometry of the doping fronts. Though the 
experiments of Ref concerned only 2D OSC films, the 
3D bulk geometry of the doping process is also techni- 
cally attainable for transparent electrodes made, e.g., of 
chemically derived grapheneiS,. In the 3D geometry, the 
stationary weakly nonlinear equation (j24p becomes 

AU/Up = (1 + AU/Up) Jf + {Vf f + ^VV, (42) 

where V-operator concerns two transverse coordinates, 
y and z. By using the methods of combustion 
theory^idSii^j we may look for a particular solution to 
Eq. in the form 
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FIG. 11; P-doping fronts with small-scale (a) and large-scale 
(b) corrugations as compared to the sample size, with periods 
A = 0.3mm and 1.2mm, respectively. Both fronts propagate 
with the same maximal velocity. 



The solution (|43p describes only rectangular 3D patterns 
at the doping fronts as illustrated in Fig. 12. Still, com- 
parison of triangular, rectangular, hexagonal and axisym- 
metric patterns for a similar combustion problem demon- 
strated that increase of the front velocity depends slightly 
on a particular patterii^. The combustion studies sug- 
gest also the evaluation max(At/3£)) « 2 max(AC/2£)), 
which would imply velocity increase by a factor of 3 for 
curved p-doping fronts in the 3D bulk OSC devices. How- 
ever, careful investigation of Eq. (|42p indicates that ve- 
locity increase of a curved p-doping front with 3D pat- 
terns is even stronger than that suggested by the above 
estimate. Indeed, by substituting the solution Eq. P5)) 
into the stationary equation ([^^ , we come to the follow- 
ing equation for the maximal velocity increase in the 3D 
geometry 



AU^D/Up 



(1 -I- AU^d/UpY 



(44) 



which looks similar to Eq. (j32l) . but, surprisingly, has 
no physical solution. The last fact implies some impor- 
tant consequences for curved doping fronts in optoelec- 
tronic devices with the 3D bulk geometry: (i) the ap- 
proach of weak nonlinearity of Eq. (|23p does not work in 
the 3D geometry, and one has to employ the fully non- 
linear equation (jl9p : (ii) presumably, the doping front 
instability is so strong in the bulk geometry, that it 
cannot be not stabilized at the nonlinear stage to any 
stationary /quasi-stationary solution. Instead of the sta- 
bilization, one may expect acceleration of the 3D humps 
at the doping fronts resembling the breakdown in light- 
ning and development of streamer a^^'^^ . Such an acceler- 
ation of the doping fronts with 3D patterns may be also 
compared to acceleration of curved combustion fronts in 
deflagration-to-detonation transitio n^^i^^ . Similar to the 
present work, combustion studies have encountered an 
example of curved fiames in channels with the propaga- 
tion regime determined by the geometry dimension: in 
open tubes, the combustion fronts propagate stationary 
in average in the 2D geometry, but accelerate in the 3D 
Thus, one should expect new interesting 



case, see 



48.49 



hoiy^z) = f2D[y) + f2D{z). 



(43) 



effects for doping front dynamics in the 3D bulk optoelec- 
tronic devices, which are qualitatively different from the 
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FIG. 12: A corrugated p-doping front with a rectangular 3D 
pattern. 

case of 2D OSC films, and which may further improve 
the device kinetics. 



VI. SUMMARY 

In the present work we have developed a nonlinear the- 
ory of the doping front instability discovered recently in 
Ref li^. Our nonlinear studies predict the possibility of 
quasi-stationary curved p-doping fronts in OSC films, 
such as the ones employed, e.g., in light-emitting cells. 
The curved doping fronts may propagate with a veloc- 
ity that is increased by a factor of 2 in comparison with 
planar fronts, thus allowing for an optimization of the 
optoelectronic devices. The theoretical predictions for 
the velocity increase are in good agreement with the ex- 
perimental measurements, see e.g. RefJ^. According to 
the theory, further increase of the front velocity may be 
achieved by fine tuning the initiation of the instability 
using minor smooth humps introduced at the electrodes. 



The theory also predicts the front shape in the form 
of smooth humps facing the undoped material and sharp 
cusps facing the doped substance, in agreement with ex- 
perimental observations. At the same time, in contrast 
to analogous studies of combustion front dynamics, ex- 
tra deep cusps at the doping fronts do not produce addi- 
tional increase of the front velocity; instead, they lead to 
certain velocity reduction in comparison with the moder- 
ately corrugated doping fronts. In the experimentally im- 
portant parameter domain, humps at the p-doping fronts 
demonstrate self-similar properties, which allows for bet- 
ter control of the front shape and the doping rate. In 
particular, we suggest that a maximal doping rate may 
be achieved both for strongly corrugated fronts, and for 
fronts which look planar on average on large scales but 
demonstrate corrugations on small scales. This interest- 
ing property of the doping front instability allows control 
of the shape of the light-emitting p-n junction. 

Finally, the present analysis suggests that dynamics 
of the curved p-doping fronts in optoelectronic devices 
with a 3D bulk geometry may be qualitatively different 
from that in the 2D OSC films studied in this paper. We 
expect that p-doping fronts with 3D patterns will not 
be stabilized to a stationary (or quasi-stationary) non- 
linear pattern even in the case of a single front propa- 
gating in an initially uniform external electric field, but 
they will rather accelerate, similar to breakdown in light- 
ning. Such an acceleration may reduce the turn-on time 
in the 3D optoelectronic devices considerably in compar- 
ison with the light-emitting cells based on 2D OSC films. 
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